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A small quantum absorption refrigerator, consisting of three qubits, is discussed in the transient 
regime. We discuss time scales for coherent dynamics, damping, and approach to the steady state, 
and we study cooling and entanglement. We observe that cooling can be enhanced in the transient 
regime, in the sense that lower temperatures can be achieved compared to the steady-state regime. 
This is a consequence of coherent dynamics, but can occur even when this dynamics is strongly 
damped by the dissipative thermal environment, and we note that precise control over couplings or 
timing is not needed to achieve enhanced cooling. We also show that the amount of entanglement 
present in the refrigerator can be much larger in the transient regime compared to the steady-state. 
These results are of relevance to future implementations of quantum thermal machines. 


Recently, the study of small self-contained quantum 
thermal machines has received growing interest; see [1- 
4] for recent reviews. Such machines typically consist 
of only a few quantum levels, hence can be considered 
'small' quantum systems (in terms of Hilbert space di¬ 
mension). Moreover, these machines are termed self- 
contained (or autonomous) as they function without any 
source of work or external control, but use only heat 
baths at different temperatures. The simplicity of these 
models makes them an ideal testbed for investigating 
quantum thermodynamics [3]. 

First works in this area go back to the study of the 
thermodynamics features of lasers [5]. Since then, many 
designs have been proposed and studied (see e.g. [6-9]), 
among these a quantum absorption refrigerator consist¬ 
ing of three qubits [10]. The efficiency of this machine 
[11] and more general performance bounds [12,13] were 
discussed. The basic functioning and fundamental lim¬ 
its of the fridge can be captured via the concept of vir¬ 
tual qubits [14]. Moreover, quantum entanglement was 
shown to appear in this model, and to enhance cooling 
in certain regimes [15]. Possibilities for experimental 
implementations [16-18] were also discussed. 

So far, most works have discussed quantum absorp¬ 
tion refrigerators in the steady-state regime, giving a 
detailed characterization of its physical properties. On 
the other hand, the transient regime remains basically 
unexplored so far. The latter is however of interest. First, 
from a conceptual point of view, it is relevant to under¬ 
stand the approach to equilibrium. Second, from a more 
applied point of view, it is natural to ask how fast cool¬ 
ing can be achieved, and what the timescale for reaching 
equilibrium is. The study of quantum effects, such as 
entanglement and coherence in the transient regime is 
also an interesting issue. 

Here we investigate the physics of a quantum absorp¬ 
tion refrigerator in the transient regime. We focus on the 
three-qubit quantum fridge model of Ref. [10], character¬ 
izing time scales, cooling properties and entanglement. 


First, the time scales for coherent dynamics, damping, 
and decay to the steady state in terms of the bath coup¬ 
ling and interaction strengths are discussed. Then we 
observe that cooling can be enhanced in the transient 
regime via the coherent dynamics of the system. Spe¬ 
cifically, it is possible to bring the object to be cooled to 
a temperature which is much lower than its steady-state 
temperature. As discussed recently in Ref. [19], this 
is a genuinely quantum feature, which highlights the 
advantage offered by quantum refrigerators over purely 
classical ones. Moreover, we find that neither precise 
timing nor accurate control of the coupling strengths 
are necessary for taking advantage of this cooling en¬ 
hancement. Finally, we also observe that the amount 
of entanglement that can be achieved in the model can 
be much larger in the transient regime compared to the 
steady-state regime, which is again a consequence of the 
coherent dynamics. We believe that the present results 
opens novel questions for quantum thermal machines 
and may be of relevance to future possible practical 
implementation of these ideas. 

I. MODEL AND MASTER EQUATION 

We consider the model of a three-qubit quantum 
absorption refrigerator discussed in Refs [io, 11] and 
sketched in Fig. 1. The ground and excited state of qubit 
i are denoted by |0) z - and |1) • and the energy gap E*. The 
free Hamiltonian for the system is thus given by 

H 0 = £ Ei |1),-(1|. (1) 

ie{C,R,H} 

We set Er = Eq + Eh and Ec 7^ Eh, and thus the states 
1010) and 1101) are degenerate in energy (we use the 
ordering CRH). Moreover we consider an interaction 
Hamiltonian 

Hint — &(| 010 ) (1011 + 1101 ) ( 010 |), 


(2) 
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Figure 1. (Color online) Three-qubit quantum thermal machine. 
Three two-level quantum systems (qubits) with energy gaps 
Ecf E r , and Eh are coupled to separate thermal reservoirs at 
temperatures Tq, Tr, and T h with coupling rates pc, PR / and 
Ph- The qubits interact collectively as described in the text, 
with interaction strength g. When operating as a fridge, the 
machines cools qubit C, i.e. brings it to a temperature below 
To 

where g is the coupling strength. We focus on the weak 
coupling regime, g <C E/, hence all transitions (apart 
from 1010) 1101)) are exponentially suppressed. 

Each qubit is weakly coupled to a thermal bath. The 
first qubit is connected to the coldest bath, at temper¬ 
ature To The third qubit is connected to the hottest 
bath, at temperature Th . The middle qubit is con¬ 
nected to a bath at intermediate temperature Tr, with 
Tc < Tr < Th- The interaction between each qubit and 
its bath is modeled by a simple reset model (see e.g. [n]) 
where thermalisation happens through rare but strong 
events. At every time step, each qubit i is either reset to 
a thermal state T\ at the temperature of its bath with a 
small probability, or left unchanged. The evolution of 
the qubit states is thus given by the master equation 

|U-i[H 0 + fW]+ E Pi(*i(p)-P) ( 3 ) 

° l ie{C,R,H} 

where is the thermalisation rate for qubit i and 

®i(p) *= Ti®Tii(p) ( 4 ) 

where Tr z denotes the partial trace over qubit i and the 
tensor product is to be taken at position i. The thermal 
states are given by T; = r z 10) (0| + (1 — r z ) |1) (1| with 


where T) is the reservoir temperature for qubit i 
(throughout the paper we set /c# = 1 ). We note that 
the master equation applies in the perturbative regime 
where pi,g <C Ec, Eh and Pi <C 1 (see e.g.[n]). In this 


case, thermalization events between a given qubit and 
the heat bath associated with the other two qubits are 
second order events which can be safely neglected. 

In [10] it was shown in detail how the machine can 
operate as a refrigerator, considering the steady-state re¬ 
gime. Briefly, when the reduced states of each qubit are 
diagonal, we can associate temperatures T c , T r , Tj 2 with 
them via (5). This is the case in the steady state, where 
the cooling effect can be simply understood in the vir¬ 
tual qubit picture developed in [14]. The interaction Hj nt 
effectively puts the cold qubit in thermal contact with 
a virtual qubit spanned by the levels |01)r H and |10)^ H 
of the other two qubits. In the absence of interaction, 
the temperature of this virtual qubit is 

Tv = ErTTr^EhTTh' (6) 

When 0 < Ty < Tq, the machine acts as a refrigerator, 
in the sense that the cold qubit will be cooled below its 
bath temperature, i.e. Ty < T c < Tq. 

II. SOLVING THE MASTER EQUATION 

Previous works have discussed the steady state 
solution of the above model in great detail, see e.g. 
[io, n, 15]. Here our focus is different as we are in¬ 
terested in the transient regime. 

We start by pointing out that (3) is linear in p and can 
thus be recast as a matrix differential equation 

9 v 

¥ = Av + u, ( 7 ) 

where v is simply a rewrapping of the density mat¬ 
rix p to a vector. The matrix A and vector u depend 
on the parameters E z , g, p- u T), and encode the right- 
hand side of (3). The steady-state solution is given by 
Voo = — A _1 u. At intermediate times, (7) is solved by 
v(f) = v h (t) + Voo, where v/ z (f) is a general solution to 
the homogeneous equation with u = 0, which can be 
obtained by diagonalising A. Denoting the eigenvalues 
and eigenvectors of A by Ay and ey respectively, one has 

v /<(0 = L c i e¥e i' ( 8 ) 

i 

where the coefficients cy are determined by the initial 
condition that v( 0 ) matches the given input state. Note 
that the real part of all A z must be negative such that 
v/ z (f) vanishes at long times and the steady state Voo is 
recovered. 

It turns out that diagonalizing A in full generality 
is challenging. Nevertheless, for fixed values of the 
parameters E/, g, pi, T), we can easily obtain the time- 
dependent solution to (7) and hence the state p(t) for any 
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given initial condition p(t = 0 ) = po- We also know the 
steady-state poo for given parameters (using the method 
above, or from [11]). 

III. CHARACTERIZING THE TRANSIENT REGIME 

We are now in position to start discussing the phys¬ 
ical properties of the three-qubit fridge in the transient 
regime. We will first consider the time scales involved 
in the approach to the steady state, and then look at 
cooling and entanglement in the transient regime. 

Throughout the following we will consider a fixed, 
natural initial state, namely a thermal state where each 
qubit is at equilibrium with its bath 

T = T C (8) Tr <g> T H . (q) 

This is the equilibrium state corresponding to the free 
Hamiltonian Hq that one has before the interaction is 
turned on. Initially there is thus no entanglement, and 
each qubit is at the temperature of its respective bath. 
Since the state is diagonal in the basis of Hq the only 
elements of the density affected by the unitary evolution, 
once the interaction is turned on, are those in the degen¬ 
erate subspace affected by Hj nt . In particular, since the 
interaction with the baths do not generate coherences 
within each individual qubit, the only off-diagonal ele¬ 
ments will be | 010 ) (1011 and 1101 ) ( 0101 . This also im¬ 
plies that the reduced state of each qubit is diagonal at 
any time, and we can thus associate temperatures with 
them. 

In addition to fixing the initial state, we will focus on 
the role played by the relative strengths of the couplings 
pi and the interaction g. Without loss of generality, we 
can set the energy scale such that Eq = 1. We also fix the 
other qubit energy and the bath temperatures. For sim¬ 
plicity we set Tq = Tr (this also analogous to a natural 
situation where the cold bath is at room temperature). 
As long as the qubits are sufficiently far from each other 
we can model their baths as independent [19, 20]. 

A. Quantities of interest 

Before moving on, we define the quantities of interest 
we will use below to study the approach to the steady 
state, the cooling, and the entanglement in the transient 
regime. 

To examine the approach to equilibrium, one needs a 
natural notion of distance from the steady state, and we 
take the trace distance 

1 

D(p(t),poo) = -|0°°||l, (10) 

which measures the distinguishability between the two 
states (|| • ||i is the trace norm). Specifically, this quantity 


has an operational meaning, as it is equal to the classical 
trace distance between the probability distribution of 
the measurement outcomes for the optimal quantum 
measurement distinguishing between p(t) and poo [21]. 

For cooling, we note that starting from the product 
thermal state r, the only off-diagonal elements of p(t) at 
any time are those in the degenerate subspace. Therefore 
the reduced state of each qubit is diagonal, and using 
(5) we can associate temperatures T C/ T r , T \ with them. 
When operating as a refrigerator the machine cools the 
cold qubit such that T c < Tq- 

Finally, for the entanglement in the three-qubit system, 
there are several bipartitions and quantities one may 
study. Entanglement in the steady state was discussed 
in Ref. [15], and all types of three-qubit entanglement 
were shown to occur in various regimes. Following this 
work, we evaluate entanglement along a given bipar¬ 
tition (e.g. qubit 1 versus qubits 2 and 3), or genuine 
tripartite entanglement, using a class of entanglement 
witnesses developed in [22, 23] which allow one to fully 
characterize the entanglement of states of the form p(t) 
for initial state r. Specifically we consider witnesses of 
the form 

w s (p(t)) = 2 ^ 1 p3,6 1 - E ^ °> (“) 

where p^ denotes elements of the density matrix p(t) 
and the set S depends on the partition and type of 
entanglement one is interested in. When inequality (11) 
is violated, its left-hand side gives the concurrence [24] 
of C|RH (S = {2}), R|CH (S = { 1 }), CR\H (S = { 3 }) 
or the genuine multipartite concurrence (see Refs. [25, 
26]) for S = { 1 , 2 , 3 }. When inequality (11) holds, no 
entanglement is present on the given bipartition, as the 
witness provides a necessary and sufficient condition for 
biseparability [27]. 

B. Time scales 

In Fig. 2 we show representative plots illustrating the 
transient behaviour of the distance to the steady state, 
the entanglement, and cooling. We observe the follow¬ 
ing general behaviour. When the interaction strength g 
exceeds the bath couplings p\, all of the observed quant¬ 
ities initially oscillate at a frequency of approximately 
g/ re (see Fig. 2). This is intuitive since this is the time 
scale of the system dynamics in the absence of dissipa¬ 
tion (as can be seen from the interaction Hamiltonian). 
Hence for weak bath coupling, we expect to see such 
oscillations until dissipation becomes dominant. 

More specifically, from (8) the time scales are determ¬ 
ined by the eigenvalues Ay of the matrix A. Since starting 
from the initial state r, p(t) always features a single off- 
diagonal element, we can take A to be 9 x 9 . We observe 
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Figure 2. (Color online) Plots vs. time of (a)-(b) the distance to the steady state D(p(t),p 0 o), (c)-(d) the cold qubit temperature T c (t), 
(e)-(f) bipartite entanglement and (g)-(h) genuine tripartite entanglement WcRH(t)- We use ^ wo se * s parameters 

given by E c = 1, E H = 100 , T c = T R = 1 , T H = 100 , pc = Vh — 10 ~ 5 , Pr = 10 -3 , and g = 10 -2 (blue curves), g = 10 -4 (red 
curves), and the system is initially in a thermal state with each qubit equilibrated to its bath. In (c) and (d) the minimal cold 
qubit temperature attainable by unitary dynamics alone is indicated (dashed horizontal), and in (d) also the period of the unitary 
dynamics (dashed vertical). In (e)-(h) the maximal entanglement extractable from the initial state by the unitary dynamics alone 
is indicated (dashed). 


the following general properties of the spectrum of A. 
Only two of the eigenvalues are complex (conjugates 
of each other) \ cp , A*^. Moreover, all eigenvalues have 
negative real parts (ensuring convergence to the steady 
state). The largest (numerically smallest) eigenvalue 
A m ax is always real. 

The coherent dynamics is characterized by oscillations, 
the frequency of which is given by the imaginary part 


of A C p, while the time scale for damping is given by the 
real part. When g is small compared to the couplings 
Pi, the real part of \ cp is pc + PR + Vh> while for large g 
it is 3 (pc + Pr + Ph)/ 4 - A typical example is shown in 
Fig. 3(a). The oscillations are therefore damped out in 
a time that scales as the inverse of pc + Pr + Ph- Thus, 
this is the time it takes for the dissipative processes to 
suppress coherent dynamics in the system (analogous 
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to a T2-time for the machine in the language of spin 
relaxation; the time to equilibrate to the steady state 
would be a Ti-time). If the sum of the bath couplings 
exceeds g, no oscillations are observed, as apparent from 
the red curves in Fig. 2, while for large g the imaginary 
part of A C p is 2 g as expected. 

After coherent dynamics is damped out, each quant¬ 
ity approaches its steady-state value. The trace distance 
initially drops fast, in the coherent regime, and then fol¬ 
lows an exponential decay. The rate of this approach to 
the steady state at long times is given by A mflX . Denoting 
the rewrapping of v to a density matrix p (v) this can be 
seen easily 

D(p(t),poo) = l||p(vfc(0+v oo )-Kvo O )||i 

= ^ IKEA'S) Hi 

i 

~ 2 I \p( c max£^ maxt £max) | |l 

= 2 Cmax ^ mait I \p( G max) | |l/ (l 2 ) 

where c max and e max are the coefficient and eigenvector 
corresponding to A mflX , and we used that the wrap¬ 
ping is linear and (8). In general, A mflx may depend 
on all the bath couplings, the interaction strength, and 
the temperatures. However, we observe that if the 
two smaller p z are equal, then \ max is equal to this 
value, X max = min {pcsPRsPh}, independent of g and 
the temperatures. This is also the value of A mflX if g 
is small compared to the p z . For large g, we have 
A max ~ Vmin + v'mJ ^ where p min , p' min denote the two 
smallest couplings. Typical examples illustrating this 
behaviour are shown in Fig. 3(b). 

We note that it can be helpful to think of the machine 
as analogous to a damped oscillator, with the regimes 
where coherent evolution is visible or suppressed cor¬ 
responding to under- and overdamping respectively. 


C. Cooling 

Both steady-state and transient cooling increases with 
increasing interaction strength. When the goal is to 
reach a low temperature for the cold qubit, it is always 
optimal to make g as large as possible. The time it 
takes to cool, however, depends strongly on whether 
the system displays coherent dynamics or not, i.e. on 
whether g exceeds the p\, and depending on the bath 
couplings, it can happen that colder temperatures are 
reached in the transient than in the steady-state regime, 
as was also noted in [19]. 

The difference between evolution with and without 
the coherent dynamics damped out is apparent from 
the examples in Fig. 2. In Fig. 4(a) we show the coldest 
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Figure 3. (Color online) (a) Damping rate of coherent dynamics 
vs. interaction strength for Eq = 1 , E# — 100/ Tc = 1, Tr = 1, 
Th = 100, pc = 10 -4 , Pr = 10 -3 , pH = 10 -4 . The limiting 
values pc + Pr + Pc an d 3 (pc + Pr + Pc)/^ f° r small and 
large g respectively are indicated (dashed), (b) Asymptotic 
decay rate vs. interaction strength for Ec = 1 , Eh — 100 , 
Tc = 1, Tr = 1, T H = 100, pc = 10- 4 , PR = 10-3, PH = 10-5 
(blue) and pc = Ph = 2 x 10 -5 (red). The limiting values 
min{ pc, Pr, Ph} = Ph an d ~ Ph + Pc/4 for the case of all p z - 
different are indicated (dashed). 


temperature of the cold qubit reached within a fixed 
evolution time as a function of g. We see that as g sur¬ 
passes the bath couplings, there is a sharp drop in the 
cold qubit temperature attained at short times. The sys¬ 
tem cools faster when the interaction is strong such that 
the initial dynamics are coherent, as might be expected. 

To separate out the effect of the dissipative evolu¬ 
tion, we can compute the lowest cold qubit temperature 
which can be reached from the initial thermal state r 
under unitary evolution alone with the Hamiltonian 
Hq + Hi n t. The minimal temperature is reached when 
the ground state population of the cold qubit is maximal. 
This population (referring to the computational basis) is 
ppi + p2,2 + P3,3 + P4,4, arid since only p 3^ is affected by 
the evolution, that means p 3^ should be maximal. The 
evolution shuffles population between the states |010) 
and 1101), corresponding to p 3^ and p6,6 and hence the 
maximal value of p3 r 3 is max {T3 r 3, Which of these 
elements is larger depends on the bath temperatures and 
the energies. However, in a regime where the machine 
provides steady-state cooling, the larger one is always 
T6,6. This can be seen using the virtual qubit picture [14]. 
Expressed in terms of the cold and (normalised) virtual 
qubit ground state populations at time zero, we have 
t 3,3 = rc(l — ry) and = (1 — rc)ry. Steady-state 


















6 



9 9 




Figure 4. (Color online) (a) Cold qubit temperature vs. in¬ 
teraction strength. The lowest cold qubit temperature at¬ 
tained within a fixed time of t = 500 is shown (solid) as 
well as the steady-state value (dashed), for Eq = 1, E # = 100, 
Tc = Tr = 1, T h = 100, p c = p H = 10“ 5 , and p R = 10“ 3 . 
(c) Cold qubit temperature vs. time corresponding to the bath 
couplings of (a) and g = 5 x 10 -3 . The lowest temperature 
attainable from the initial state by coherent dynamics is indic¬ 
ated (dotted) as well as the steady-state value (dashed), (b) 
and (d) Same as (a) and (c), with pc = 10 -4 . Interestingly, the 
transient temperature of the cold qubit is significantly lower 
compared to the steady-state value. 


cooling requires Ty < Tq [ 14 ] and hence tq < ry. Thus 
the minimal cold qubit temperature is reached when the 
populations of 1010 ) and 1101 ) are completely reversed, 
which happens after half a period at time t = rc/lg. 

Interestingly, due to the coherent dynamics, the cold 
qubit temperature in the transient regime can be signi¬ 
ficantly lower than the steady state value. As the initial 
dynamics is effectively dissipation free, the minimal 
temperature in this regime is roughly independent of 
the bath couplings, contrary to the steady-state value. 
An example of this behaviour is shown in Fig. 4 . In 
particular Fig. 4 (d) shows how a cold temperature is 
quickly reached during the coherent phase of the evol¬ 
ution, while after damping out of the oscillations, the 
temperature then approaches a higher steady-state value. 
In this case, optimal cooling is achieved at t = rc/lg. 

While this optimal cooling would require a precise 
timing, we note however that even if control on this 
short time scale is not available, it may still be possible to 
achieve enhanced cooling by extracting the qubit before 
reaching the steady state. While coherent oscillations are 
damped out at a time set by the largest bath coupling, 
the approach rate to the steady state is determined by 
the smaller bath couplings, as discussed above. Hence 



Figure 5. (Color online) A minimum cold qubit temperature 
is reached in finite time for Eq = 1, Eh = 100, Tc = Tr = 1, 
T h = 100, Pc = 10- 5 , p R = 1(T 3 , p H = 10- 5 , g = 1(T 4 . 
The steady-state temperature is indicated (dashed) and half a 
period zr/2 g of coherent evolution without dissipation (dot¬ 
ted). 


when the couplings are different, there is an interme¬ 
diate regime between coherent dynamics and steady 
state where the qubit temperature is low and precise 
time control is not necessary to extract it (c.f. Fig. 4 (d)). 
A minimum cold qubit temperature can occur in finite 
time even when coherent dynamics is not visible. This 
can happen e.g. when the coupling to the hot bath is 
very weak, while the the largest bath coupling exceeds 
the interaction strength, damping out oscillations. An 
example is shown in Fig. 5 . The optimal time in this case 
scales linearly with p^ and inversely with g but occurs 
much later than the half-period rc/2g corresponding to 
the first minimum under coherent dynamics. One can 
understand this, as well as the case where some initial 
oscillations are visible, in the picture of an underdamped 
oscillator approaching critical damping. As the damping 
is increased, in addition to the amplitude of oscillations 
decreasing their period increases and finally diverges as 
critical damping is reached. 

We note that the effect of reaching a minimal temper¬ 
ature in finite time is robust to small deviations in the 
coupling and interaction parameters. The qualitative 
behaviour persists as long as the ordering of their mag¬ 
nitudes is preserved. Hence very precise control over 
the couplings is also not required to achieve enhanced 
cooling in the transient regime. 


D. Entanglement in the transient regime 

We can thus characterize the dynamics of entangle¬ 
ment. Here we focus on genuine tripartite entanglement, 
and on entanglement on the bipartition R\CH, i.e. the 
bipartition of energy in vs energy out. Examples are 
given in Fig. 2 (e)-(h). It is clearly seen that in the coher¬ 
ent regime, the amount of both types of entanglement 
can be considerably larger than in the steady state. 

In fact, it is simple to estimate how much entan- 
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glement can be created in the transient regime. This 
amount will depend on the relative sizes of the interac¬ 
tion strength and the bath couplings, as well as on the 
initial state. However if the interaction strength is larger 
than the bath couplings, at short times the evolution 
is approximately dissipation free, as mentioned above. 
The system thus evolves approximately unitarily under 
the Hamiltonian Hq + Hj nt , and the maximal entangle¬ 
ment is what can be generated from the initial state r 
under this evolution. The only elements of r which can 
change are those in the degenerate subspace affected by 
Hj nt . The second term in (11) is then invariant and the 
maximal entanglement is obtained by maximising \ps f e\ 
which occurs after a quarter period at t = n/4g. The 
optimal value is \p3 f e\ = \r^ — r^\ which results in a 
maximal bipartite entanglement of 


ujmax _ at 
Vy/ R\CH ~ iN/ 


-li 


e c+ e h 
> t r 


E C I e h 

-e^ +T ^\ 


— 2N -1 e 2 V r c 


'Tr, 


e H ( 1 


1 


2 \T R ^T H ; 


(13) 


where N = (1 + e Ec ^ Tc )(l + c Er / Tr )(1 + e En ^ Tn ). The 

maximum is the same across the other bipartitions, how¬ 
ever we note that once the thermal dissipation is accoun¬ 
ted for, entanglement on R\CH appears to dominate the 
other two. For genuine tripartite entanglement, Wj^^ 
is obtained from (13) by replacing replacing the factor 2 
with 6 in the negative term. 

As an example, the maximal bipartite and genuine 
tripartite entanglement extractable by the unitary dy¬ 
namics is indicated on Fig. 2(e)-(h). It represents an 
upper bound on the entanglement which can be extrac¬ 
ted when one has control over the system at the fast 
time scale set by g. As can be seen, the maximal value 
is never reached exactly, due to the presence of the bath 
couplings. While entanglement in the transient regime 
is maximised for interaction strength large compared 


to the bath couplings, this is not optimal for generation 
of steady-state entanglement. Unlike for cooling, it is 
not optimal to make g as large as possible. As a result, 
when optimising for steady-state entanglement, no os¬ 
cillations are observed in the transient regime, and the 
entanglement remains small compared with the largest 
value obtainable from dissipation free evolution. 


IV. CONCLUSION 

The dynamics of a small quantum absorption refri¬ 
gerator was investigated. We discussed the approach 
to equilibrium, in particular with respect to time scales, 
cooling, and entanglement. Notably, we observed that 
optimal cooling may occur at finite times, and not in 
the steady-state regime. Similarly, the largest amount 
of entanglement is usually obtained in the transient re¬ 
gime. As pointed recently in Ref. [19], this highlights the 
relevance of quantum effects in such thermal machines. 

We believe that the observation that enhanced cool¬ 
ing can be achieved in the transient regime without the 
need for precise timing or control of the fridge para¬ 
meters may open interesting possibilities for cooling or 
initializing quantum systems. 

More generally, the methods discussed in the present 
work could be adapted to other models of quantum refri¬ 
gerators, as well as to other quantum thermal machines, 
producing work or entanglement [28]. 
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